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Abstract 

We analyze the electronic structure of molecules which may exist in gas phase of chemical 
vapor deposition process for GeSbTe alloy using the electronic stress tensor, with special focus 
on the chemical bonds between Ge, Sb and Te atoms. We find that, from the viewpoint of the 
electronic stress tensor, they have intermediate properties between alkali metals and hydrocarbon 
molecules. We also study the correlation between the bond order which is defined based on the 
electronic stress tensor, and energy-related quantities. We find that the correlation with the bond 
dissociation energy is not so strong while one with the force constant is very strong. We interpret 
these results in terms of the energy density on the “Lagrange surface”, which is considered to 
define the boundary surface of atoms in a molecule in the framework of the electronic stress tensor 
analysis. 
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I. INTRODUCTION 


GeSbTe (GST) alloy is the most popular material for phase change memory (PGM) 
, I 2 I], which is one of the most promising candidates for the next-generation memory device. 
So far, GST thin dims are deposited by physical vapor deposition such as guttering and 
pulsed laser deposition. However, chemical vapor deposition (GVD) of GST has many 
advantages such as good step coverage, uniformity and high purity, and considered to be 
necessary for future PGM applications. The GVD process for GST is relatively a new held 
of research and there remains many things which are not clearly known. One of them is the 
chemical reactions in the gas phase of GVD process, and we have investigated reactions and 
molecules which may exist in the gas phase using quantum chemical calculation jb]. As the 
result of this study, we have obtained a data set of molecules which have bonds among Ge, 
Sb and Te. In this paper, we apply our electronic stress tensor analysis based on the Rigged 
QED (Quantum Electrodynamics) theory [t-ISI to this data set, and investigate how these 
bonds can be described by the electronic stress tensor. 

In general, the stress tensor, which describes a pattern of internal forces of matter, is 
widely used in various helds of science such as mechanical engineering and material science. 
The use of stress tensor in quantum systems as well has been investigated for many years, by 
one of the earliest quantum mechanics papers IJ, [iSj and many researchers 16l-l3l| including 
jr-lS, 32-42\. In our past studies, we have shown that the electronic stress tensor 


our group 


and related quantities can be useful tools to analyze atomic and molecular systems and can 
give new viewpoints on the nature of chemical bonding. 

We here introduce two of our Endings which we wish to investigate more deeply using GST 
bonds in this paper. First, we have pointed out that we may characterize some aspects of a 
metallic bond or metallicity of a chemical bond in terms of the electronic stress tensor 32l.l42|. 
We have analyzed the bonding region in the small cluster models and periodic models of 
Li and Na using the electronic stress tensor, and have found that all the three eigenvalues 
of the stress tensor are negative and degenerate, just like those of liquid. This is in stark 
contrast to hydrocarbon molecules, which have the positive largest eigenvalue much larger 
than the other two negative eigenvalues. The eigenvalue pattern of Li and Na indicates a 
lack of directionality and compressive nature of bonding while that of hydrocarbon molecules 
indicates solid directionality and tensile nature of bonding. Each pattern well reflects the 
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nature of metallic and covalent bonding. Then, the question worth asking is, how the 
chemical bonds between metalloid atoms like Ge, Sb and Te are described by the electronic 
stress tensor. Second, using the energy density which is defined from the electronic stress 


tensor, new dehnition of bond order has been proposed 


a correlation between this new bond order and bond distance 


So 

ar, 

we 

32 

34, 

0 


that our bond order exhibits better correlation than the other bond orders proposed in 
the literature. We, however, have not investigated a correlation with a quantity related to 
energy. Therefore, we wish to investigate a correlation between our bond order and the 
bond dissociation energy, which is available in our data set as the dimerization energy of the 
CVD precursors. In addition, we compute a force constant for the GST bonds and examine 
a correlation with our bond order. 

The structure of this paper is as follows. Sec. [IT] summarizes our analysis method of 
electronic structures using the electronic stress tensor density. In Sec. m we show the 
electronic stress tensor density and its eigenvalues of the GST bonds, and compare with 
those of hydrocarbon molecules and alkali metal clusters. We also investigate correlations 
between our bond order and energy-related quantities. Finally, Sec. HV] is devoted to our 
conclusion. 

II. THEORY 


In the following section, we analyze chemical bonds using quantities such as the electronic 
stress tensor density and the kinetic energy density. They are based on the Rigged QED 
theory Sj and we briefly describe them in this section. The Rigged QED is a quantum field 
theory which has been proposed 8 | to treat dynamics of charged particles and photons in 
atomic and molecular systems. In addition to the ordinary QED which contains the Dirac 
held for electrons and the gauge held for photons, the Schrodin ger helds for atomic nuclei are 


included. More details are found in our previous papers 


7|,l8|, 


11 
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Below, c denotes the 


speed of light in vacuum, h the reduced Planck constant, e the electron charge magnitude 
(so that e is positive), and m the electron mass. The gamma matrices are denoted by 7 ^ 
(/i =0-3). 

The most basic quantity in the Rigged QED is the electronic stress tensor density operator 
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jg (defined as follows 
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jikif \ ^ 

^ 2 


^jJ{x)'J^Dek{x)^jJ{x) - ( Dek{x)^p{x)) 'y°'y^ijj{x) 


( 1 ) 


where 'ijj{x) is the four-component Dirac field operator for electrons, the dagger as a super¬ 
script is used to express Hermite conjugate, and '0(x) = 'ip\x)'y^. We denote the spacetime 
coordinate as x = {ct, r). The Latin letter indices like k and I express space coordinates from 
1 to 3. Here, the gauge covariant derivative is dehned by Dek{x) = dk + i^Ak{x), where 
Ze = —1, and Ak{x) is the vector potential of the photon held operator in the Coulomb 
gauge (divH(x) = 0). The important property of this quantity is that the time derivative of 
the electronic kinetic momentum density operator ne(x) 

Ae{,x) = ^ (^h'^\x)De{x)'^{x) - ik j , (2) 


can be expressed by the sum of the Lorentz force density operator Le(x) and the tension 
density operator ff^(x), which is the divergence of the stress tensor density operator: 


Wt 


YiJx) = Leix) + fe 


X 


( 3 ) 


These operators are expressed as follows. 


Leix) = Eix)peix) -jeix) X H(x), 

c 


rf\x) = difr\x) 


Tlkl, 


( 4 ) 

( 5 ) 


ihc 
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{belix)-li)ix)^ 


qV ■ -Defc(a:)^/’(x) -F i>ix)'y^Dekix)Deiix)i>ix) 


- (^L)ekix)L)elix)i!ix)y ■ '0(x) - (a^)) S V ■ De;(x)V^(x) 

(jeix) X H(x)) , (6) 


where i?(x) and -B(x) denote the electric held operator and magnetic held operator respec¬ 
tively, and Peix) and jeix) are the electronic charge density operator and charge current 
density operator respectively. 
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For nonrelativistic systems, we approximate the expressions above in the framework 
of the Primary Rigged QED approximation 1^ 131, in which the small components of 
the fonr-component electron held are expressed by the large components as ~ 

— ^^iha^DkilJLix) and the spin-dependent terms are ignored. Then, we take the expec¬ 
tation valne of Eq. ([3]) with respect to the stationary state of the electrostatic Hamiltonian. 
This leads to the eqnilibrinm eqnation as 


0 = (Te(a;)) + (r^ix)) = {L^{x)) + di{f^'^\x)), 


(7) 


which shows the balance between electromagnetic force and qnantnm held force at each 
point in space. Since this expresses the fact that the latter force keeps the electrons in the 
stationary bonnd state in atomic and molecnlar systems, we can expect to stndy the natnre 
of chemical bonding from the viewpoint of qnantnm held theory by nsing the stress tensor 
density and tension density. We express (rf (x)) and {f^^\x)) respectively and 

for simplicity (we also write only spatial coordinate r becanse we consider stationary state). 
The explicit expression for the stress tensor density and tension density r'^^(r) are 



4m ^ * * dx^dx^ dx^ dx^ 

i L 

dx^dx^ * ^ dx^ dx^ 


( 8 ) 
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dxk Qxk 

dA'i/j* (f) 


dx^ 


V’i(r) - 


d'lpijr) 

dx^ 


(9) 


where is the ith natnral orbital and Vi is its occnpation nnmber. A denotes the 

Laplacian, A = Y^k=ii^/■ When the density fnnctional theory (DFT) method is nsed 

to compnte the electronic structnre, we nse the Kohn-Sham orbitals for in the above 

o-s 

expressions. The eigenvalne of the symmetric tensor r is the principal stress and the 
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eigenvector is the principal axis as follows: 
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We use a concept of “Lagrange point’' 


32l | to characterize a bond between two atoms. 
The Lagrange point is dehned as the point where the tension density r‘^(r^ vanishes, 
namely T^^{rL) = 0. We analyze chemical bonds by computing the eigenvalues of electronic 
stress tensor density at this point. 

Another important quantity in the Rigged QED is the electronic kinetic energy density 
operator dehned as [3|, 


fe{x) = - 


2m 




( 12 ) 


As is done for the electronic stress tensor density operator, we apply the Primary Rigged 
QED approximation to Eq. (IT^ and take the expectation value with respect to the stationary 
state of the electrostatic Hamiltonian. Then, we obtain the dehnition for the electronic 
kinetic energy density as 

^Te(fO = ■ (13) 


Note that our dehnition of the electronic kinetic energy density is not positive-dehnite. Using 
this kinetic energy density, we can divide the whole space into three types of region: the 
electronic drop region Rjj with > 0, where classically allowed motion of electron is 

guaranteed and the electron density is amply accumulated; the electronic atmosphere region 
Ra with nT^{f) < 0, where the motion of electron is classically forbidden and the electron 
density is dried up; and the electronic interface S with nT^{r) = 0, the boundary between 
Ro and Ra, which corresponds to a turning point. The S can give a clear image of the 
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intrinsic shape of atoms and molecules and is, therefore, an important region in particular. 

Finally, in our analysis, we use the energy density sf{r) which is dehned as a half of the 
trace of the stress tensor density [^: 



k=l 


(14) 


It can be regarded as the energy density in a sense that the integration over whole space 
gives usual total energy E of the system: J e^{r)dr = E. This can be proved by using the 
virial theorem . Using this energy density a new dehnition of the bond order (bond 

strength index) is proposed [32|. It is dehned as ef{r) at the Lagrange point between two 
atoms. Our dehnition of bond order between atoms A and B is 


^ gfAB(^L) 


(15) 


One should note we normalize by the value of a II 2 molecule calculated at the same level of 
theory (including method and basis set). 


III. RESULTS AND DISCUSSION 


A. Data set and computational details 


Our data set consists of 35 molecules (Table [H the leftmost column) which may exist in 
the gas phase of CVD process for GeSbTe alloy. They could be formed by reactions among 
GST precursors and/or II 2 carrier gas, and each of them has a bond between Ge, Sb or Te 
atoms. Their geometries and coordinates are given in the supplementary material (Fig. SI 
and Table SI). For the precursors, we assume tert-Butylgermanium (GeH 3 (tBu), Fig. [T](a)) 


for the Ge precursor, triisopropylantimony (Sb(iPr) 3 , Fig. [11(b)) 


n 


or the Sb precursor, and 


diisopropyltellurium (Te(iPr) 2 , Fig. [11(c)) for the Te precursor [Sj]. (Their coordinates are 
given in the supplementary material. Table S2.) They are optimized by the DFT method 


based on the Lee-Yang-Parr gradient-corrected functional [^ with Becke’s three hybrid 
parameters ^ (B3LYP). Threshold for maximum force is set to 0.000450 hartree/bohr. 


The Dunning-Huzinaga double-zeta basis set 46|] with effective core potential by Hay and 
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Wadt 47H49| (LanL2DZ) are used for Ge, Sb and Te atoms. 6-31G(d) 50|, 1 m| basis set is 
used for G atoms. D95(d,p) basis set is employed for H atoms. Total energies at OK is 


obtained including zero-point energies with a scalin g fa ctor for B3LYP, 0.980 


We also use Ge, Sb and Te crystal structures 


53 




55l | as a reference data (see Table S3 


m 


the supplementary material for the details). We use the norm-conserving pseudopotentials of 
Troullier-Martins typ e 56| and the generalized-gradient approximation method by Perdew- 


Burke-Ernzerhof 


57( 1 for density functional exchange-correlation interactions. Kinetic energy 


cutoff of plane-wave expansion (k-point) is taken as 40.0 hartree (2x2x2 k-point set). 


The electronic structures used in this paper are obtained by Gaussian 09 


for cluster 


models and by ABINIT 59|, [GOJ for periodic models. We use the QEDynamics package 61 | 
developed in our group to compute the quantities described in the previous section such as 
Eqs. (IH]), ((9]) and ([13]). 


B. Electronic stress tensor and its eigenvalues 

We begin by briefly reviewing our past works on the electronic stress tensor analysis 
which are related to the present paper. First, it has been proposed that a covalent bond 
can be described by the eigenvalues and eigenvectors of the electronic stress tensor [9|. In 
detail, the bonding region with covalency can be characterized and visualized by the “spindle 
structure”, where the largest eigenvalue of the electronic stress tensor is positive and the 
corresponding eigenvectors form a bundle of flow lines that connects nuclei. As an example, 
we show a map of the largest eigenvalue of the electronic stress tensor including a region 
between G atoms of GeH 3 (tBu) in Fig. [2] (a), where we can find the spindle structure. In 
passing, in Fig. |2] (b), we show the tension density and the Lagrange point for the same 
G-G bond. Then, we have proposed that the negativity of the three eigenvalues of the 
stress tensor and their degeneracy, which is the same pattern as liquid, can characterize 


some aspects of the metallic nature of chemical bonding 


m 


42]. In Ref. 42], it has been 


shown that the three eigenvalues of the Li and Na clusters have almost same values while 
the hydrocarbon molecules have the largest eigenvalue much larger than the second largest 
eigenvalue, which has similar value to the smallest eigenvalue. In terms of the differential 
eigenvalues, the Li and Na clusters have very small and which are 

much smaller than of hydrocarbons. The former degeneracy pattern indicates 
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that the bonds are not directional while the latter indicates the clear directionality of the 
bonds, reflecting the metallic nature of chemical bonding in the alkali metal clusters and 
the covalent nature of bonding in the hydrocarbon molecules. 

Now, let us move on to study the electronic stress tensor of the GST bonds. We hrst show 
a map of and corresponding eigenvector on a plane which contains a bond between GST 
atoms in FigJSl We hnd a Lagrange point between each bond and its position is marked in 
the hgure. A common feature we see in all the six panels is that the eigenvectors form a 
pattern which connects GST nuclei. However, not all of them are spindle structures. As for 
the Ge-Ge (panel (a)) and Ge-Sb (panel (d)) bonds, since they do not exhibit a positive 
region, they are called to have pseudo-spindle structures |^. As for the Te-Te (panel (c)) 
and Sb-Te (panel (f)) bonds, although we may say they have spindle structures, the positive 
regions are not as conspicuous as the spindle structure of the G-G bond seen in Fig. |2] (a). 
As for the Sb-Sb (panel (b)) and Ge-Te (panel (e)) bonds, the positive regions are 
even smaller. These results lead us to conclude that the GST atoms have the ability to 
form the spindle structure in the order of Te > Sb > Ge. From our viewpoint that the 
spindle structure is the manifestation of the covalency of chemical bonding, Te contributes 
to covalency more than Sb or Ge, but less than G. The eigenvalue and eigenvector maps 
for the other GST molecules, which are found in the supplementary material (Fig. S2), also 
support this ordering. 

In Fig. m we plot at the Lagrange point against the bond distance for all the GST 
molecules. The original data are found in Table [H We see that rf> 0 for the Te-Te bonds 
and < 0 for the Ge-Ge bonds. As for the Sb-Sb bonds, can be both positive and 
negative, and absolute values are smaller than those of the Ge-Ge and Te-Te bonds, 
of the Ge-Sb bonds exhibit intermediate values between the Ge-Ge and Sb-Sb bonds, and 
similarly for the Ge-Te and Sb-Te bonds. Therefore, Fig. H] can be interpreted that the GST 
atoms contribute to the positivity of the at the Lagrange point in the order of Te > 
Sb > Ge, which is consistent with the tendency to form the spindle structure as mentioned 
above. 

We next examine the differential eigenvalues of the electronic stress tensor at the Lagrange 
points. In Fig. |5] (a), we show a scatter plot of and for GST bonds, 

and, in Fig. |5] (b), we in addition plot points for the hydrocarbon molecules, Li clusters, 
and Na clusters, which are taken from the data set studied in our previous paper |^. We 
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see that both and of the Li and Na clnsters are 0(10“'^). As for 

the GST molecules, is 0(10“^), which is somewhat larger than 

which is 0(10“^)-0(10“^). As for the hydrocarbon molecules, is 0(10“^) and 

this is much larger than Thus, the degree of degeneracy can be summarized 

as Li, Na <C GST <C h/c. This is consistent with the usual classihcation of Ge, Sb and Te 
as metalloids, which are placed in between metals and non-metals. Further research may 
reveal that the electronic stress tensor density provides a new criterion to dehne metalloids 
based on the electronic structures. 

Incidentally, we compute the electronic stress tensor of Ge, Sb, and Te crystals using 
periodic models. We compute the electronic stress tensor at the midpoint of two nearest 
neighborhood atoms (it is the Lagrange point by symmetry). The results are summarized in 
Table inland plotted in Figs. HI and [5] (a). We see in Fig. [5] (a) that all the crystal structures 
have similar values of and to those of the GST molecules. This 

indicates that the degree of degeneracy does not differ much between the crystal structures 
and molecules. As for as shown in Fig. 01 all the crystal structures exhibit negative 
values of about —2 x 10“^. This is close to the of the Ge-Ge bonds in the molecules, but 
not to those of the Sb-Sb or Te-Te bonds, which, respectively, are negative with absolute 
values of 0(10“®) or positive with values of 0(10“®). Thus, from the viewpoint of the 
electronic stress tensor, the covalency of the chemical bonding involving Sb and Te atoms 
in the GST molecules does not appear in their crystal structures and some metallicity is 
manifested. In other words, while Ge, Sb and Te do not exhibit covalency in their crystal 
structure, Sb and Te show some covalency in their molecular structure. 


C. Bond order and force constant 


We briefly review our past works 3214841. l40| on our bond order be (Eq. (IT5ll), w hich is 


dehned using concepts based on the Rigged QED. First, it has been pointed out 


32 | that be 


of a single, double, and triple bond between carbon atoms in hydrocarbons is close to 1, 2 
and 3 respectively, consistent with a conventional bond order. Also, be of the G-G bond in a 
benzene molecule is close to 1.5. However, it has been also reported that be of some diatomic 
molecules overestimate or underestimate the conventional bond order, e.g., be of N 2 is 7.462. 
Then, in Ref. 33|], the correlation between be and bond distance has been investigated 


10 









using various simple organic compounds. By comparing with other bond orders proposed 
in the literature, be is found to have a comparable or better correlation. This has been also 


shown using hydrogenated Pt [3^ and Pd 


40| clusters. We, however, have not investigated 


a correlation between and a quantity related to energy. Therefore, we here compute the 
bond dissociation energy B and force constant k of the GST bonds, and investigate the 
correlations with b^. 

The results are summarized in Table HTTl k is calculated based on the directional derivative 
of a total energy of a molecule with respect to the bond direction s. Namely, it is computed 


using force constant matrix 


d'^E 

dx'- dxi 


as 


k = 


d‘^E 

dsds 


3 


d'^E 


dx^dx^ 


( 16 ) 


The scatter plots of b^ versus r^, B, and k are shown in Fig. El We can see that b^ is negatively 
correlated with and positively correlated with k and S, consistently with the usual notion 
of bond order. In fact, their correlation coefficients are —0.968, 0.570, and 0.910, respectively 
for Te, B, and k. Since we need to know the energy when two fragments of a molecule are 
inhnitely apart to compute B, it is reasonable not to hnd a strong correlation between b^ 
and B. Note that be is dehned only using the quantities at the equilibrium structure. On 
the other hand, the reason why the correlation with k is very strong may be understood as 
follows. 

The meaning of A;, often called a spring constant, is how much energy we need to stretch 
the bond by an inhnitesimal distance. It can be expressed as k5‘^/2, where 6 denotes the 
inhnitesimal displacement. We may interpret this energy using our energy density s^{r) 
(Eq. (ITTjl ) and the “Lagrange surface” UJ, ll2| which is dehned as a separatrix in the tension 
density r‘^(r) (Eq. (I9l)). The vector held of r'^(r) generally has a pattern in which vectors 
originate from atomic nuclei, and they collide to form separatrix. See Figs. [2] (b) or [7] (a). 
We call this separatrix the Lagrange surface and regard it as the boundary of atoms in a 
molecule. In Fig. [7] (a), we show some examples of the Lagrange surface in a GeIl 3 -SbH 2 
molecule. When we move apart two atoms bounded by the Lagrange surface for a small 
distance, it may be reasonable to suppose the required energy to be proportional to the 
energy stored in the Lagrange surface, that is, where the integration is taken 

over the Lagrange surface S. To support this idea, we compute an alternative bond order 
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39l. ko ] defined as 


be{S) 




(17) 


where 5 ab denotes the Lagrange surface between atoms A and B, and investigate the cor¬ 
relation with k (Fig. S3 (c)). The correlation coefficient is found to be 0.887, and the linear 
fit y = ax + b, where x and y, respectively, represent k and &e(s), gives a = 5.54 and 
b = 5.62 X 10“^. (The value of for each bond is summarized in Table HTTl i This suggests 
that k oc e^((f)d^a holds to a large extent. As for the linear £t for b^ against k, we obtain 
a = 1.63 and b = —7.54 x 10“^. Since b is close to zero for both cases, we can approximately 
consider bs{s) oc be with a proportional constant of 3.4 (the ratio of a’s) for these GST bonds. 
This relation between be(s) and be can be confirmed directly by computing be(s)/be for each 
GST bond which is shown in Fig. [HI The average of be(s)/be over 35 bonds is 3.74 with the 
standard deviation of 0.21. Since the standard deviation is relatively small (5.6% of the 
average), we may regard be(s) oc be for these bonds. The values of the proportional constant 
derived in two ways are consistent. 

This relation, that be{s) of the GST bonds is approximately obtained by multiplying be 
by a common factor, holds if the integration of sf{r) over their Lagrange surface is well 
approximated by e^ifL) multiplied by a common factor. Such an assumption can be valid, 
if the Lagrange surfaces of the GST bonds are flat and the energy density distributions on 
them are expressed by Gaussian functions with a common value of the exponent. We can see 
this is roughly true as follows. First, the flatness can be checked by the visual inspection of 
the Lagrange surface (Fig. [7|(a) and Fig. S2 in the supplementary material). Then, we plot 
s^{f) against the distance from for the points which constitute the Lagrange surface. The 
example of this plot for the GeH 3 -SbH 2 molecule is found in Fig. [7| (b). Actually, we plot 
e^{f) divided by the energy density of a hydrogen molecule at its Lagrange point so that the 
value at tl [i.e. \f— FlI = 0) becomes be- In the hgure, we also plot the result of the £t to 
the Gaussian function fegcxpl—a|F— Fip} where the exponent a is the htting parameter. 
The figure shows that (r )/is well htted by the Gaussian function with a = 1.32. 
We perform such a £t to all the GST bonds (Figures similar to Fig. [7| (b) are plotted in 
Fig. S4. The values of a are summarized in Table IIIII and Fig. S5) and the average and 
standard deviation of the exponent are computed to be 1.31 and 9.40 x 10“^, respectively. 


12 






Therefore, as for the GST bonds, we can say that the energy density distribution over the 
Lagrange surface is expressed by the Gaussian functions with a common exponent, which 
leads to J^ef{a)(Pa(xef{fL). 

As a further test for the idea that b^(^s) is proportional to the force constant, we examine 
homonuclear diatomic molecules, from H 2 to G excluding the group 18 elements. We use 
experimental values for Vg, B, and k as is summarized in Table IIVI B is the energy of the 
ground state atomic products relative to the lowest existing level of the molecule, often 
denoted hj Dq. k is converted from the frequency u (which is also shown in the table) 
via k = Ma;^/2, where M is the atomic mass. Note that some of the experimental values 
are not available or controversial and we use computational values in such cases. We adopt 


computational values for and oj of Ga 2 from Ref. 6 ^, and and u of In 2 from Ref. [^. bg 
and bg(^s) are computed with the same setups as the GST bonds, and shown in Table HVl We 
hrst study how much b^ is correlated with r^, B, and k. The correlation coefficients are found 
to be —0.583, 0.834, and 0.978, respectively for Vg, B, and k, showing very strong correlation 
between bg and k as is the case with the GST bonds. As for 5e(s), the correlation with k is 
also very strong, with the correlation coefficient 0.990. The scatter plots of k versus bg and 
bg(^s) are shown in Fig [9] (a) and (b) respectively, and we can see these strong correlations 
(similar scatter plots for Vg and B are shown in Figs. S 6 and S7). However, one may notice 
a difference between Fig [9] (a) and (b) that the panel (b) shows more proportionality than 
the panel (a). In fact, while the linear fit bg = aik + bi gives oi = 3.12 and bi = —0.23, 
be{s) = 0 ' 2 k + &2 gives 02 = 3.64 and &2 = 5.87 x 10“^. Since b\ is not negligible, bg is hardly 
considered to be proportional to k, whereas very small value of 62 implies that bg(^s) oc k 
holds to a large extent. Therefore, in general, bgi^s) is a better descriptor of a chemical bond 
than bg. The drawback is that, since the computation of bg(^s) involves the integration over 
the Lagrange surface, it costs much greater than that of bg. 

Finally, some comments on the correlations of Vg and B with our bond orders are in order. 
First, the correlation coefficient between Vg and bg is —0.968 for the GST molecules and 
—0.583 for the homonuclear diatomic molecules. (As for 5 e( 5 ), they are —0.970 and —0.651, 
respectively.) The difference can be attributed to the fact that the GST data consists 
of similar chemical bonds formed by Ge, Sb, and Te, whereas the homonuclear diatomic 
molecule data contains various types of elements. It has been found that bg correlates well 
with Tg, as mentioned in the beginning of this subsection, but it has been also shown that the 
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slope on the re-feg plane depends on the type of chemical bonding. For example, in Ref. 40| . 
the Pd-Pd bonds and Pd-H bonds are found to have different slopes. Also, among the Pd-H 
bonds, the bonds with shorter bond length (< 1.9 A) have steeper slope than those with 
longer bond length. Therefore, when a data contains different types of chemical bonding, 
even though re and 6^ exhibit negative correlation, the correlation coefficient tends to be 
somewhat away from —1. Second, the correlation coefficient between B and is 0.570 for 
the GST data and 0.834 for the homonuclear diatomic molecules data. (As for 5£(5), they 
are 0.569 and 0.906 respectively.) There is again some notable difference between the two 
data sets. Since the correlation coefficients between and k are close to 1 for both data 
sets, the degrees of correlation between B and b^ are roughly equal to those between B and 
k. Although B and k are likely to be positively correlated, since they are not linear to each 
other in general, we do not expect a very strong correlation between B and k. Then, the 
relatively low value of correlation coefficient of the GST data is reasonable while that of the 
homonuclear diatomic molecules data is unexpectedly high. The latter may be understood 
by a universality of potential energy curve for diatomic molecules (e.g. 75l-l77l| ). which has 
been studied for a long time. It is, however, an empirical relation at this stage and the 
discussion of its relevance here is beyond the scope of this paper. 


IV. CONCLUSIONS 

We have analyzed the electronic structure of 35 molecules which may exist in gas phase 
of GVD process for GeSbTe alloy using the electronic stress tensor, with special focus on the 
chemical bonds between Ge, Sb and Te atoms. Our study consists of two parts. First, we 
have studied the pattern of the eigenvalues and eigenvectors of the electronic stress tensor 
density of the GST bonds. Next, we have computed the bond order which is defined by the 
stress-tensor-based energy density for the GST bonds, and have investigated the correlation 
with the energy-related quantities such as the bond dissociation energy and force constant. 

In the first part, we have found that, from the viewpoint of the electronic stress tensor 
density, GST bonds exhibit intermediate properties between alkali metals and hydrocarbon 
molecules. This is illustrated by the sign and degeneracy pattern of the three eigenvalues of 
the electronic stress tensor density at the Lagrange point between two atoms. In our previous 
studies, we have pointed out that the negative and degenerate eigenvalues, which indicate 
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a lack of directionality with the compressive stress, characterize some aspects of metallicity 
of chemical bonding. By contrast, the positive largest eigenvalue and much smaller negative 
other two eigenvalues, which indicates a solid directionality with the tensile stress in that 
direction, characterize covalency. The former pattern has been typically found in the alkali 
metals and the latter in the hydrocarbon molecules. Our results in the present paper suggest 
that the GST bonds can be located in between metallic bonding and covalent bonding in 
terms of their electronic stress tensor. This is consistent with the usual classihcation of Ge, 
Sb and Te as metalloids, which have intermediate properties between metals and nonmetals. 

In the second part, we have found that the correlation of our bond order with the bond 
dissociation energy is not so strong, while one with the force constant is very strong. We 
have interpreted this results in terms of the energy density on the “Lagrange surface”, which 
is considered to dehne the boundary surface of atoms in a molecule in the framework of the 
electronic stress tensor analysis. In this study, we have found that both of our definitions 
of bond order be and &e( 5 ), where the former uses the energy density at the Lagrange point 
while the latter involves the integration over the Lagrange surface, have strong correlation 
with the force constant. We have argued that, if the interpretation above is correct, he(s) is 
the one which is more directly connected with the force constant, and the strong correlation 
of the force constant with he follows from that with he{s)- In fact, we have shown that be{s)/be 
does not vary much among the GST bonds, which originate from the fact that the energy 
distributions on the Lagrange surface of the GST bonds can be well expressed by Gaussian 
functions centered at the Lagrange point and with a common value of the exponent. As the 
results of this study, it is hinted that the stress-tensor-based energy density can be related 
not only to the total energy but also to the force constant by combining with another Rigged 
QED concept, the Lagrange surface. 

In our future work, regarding the hrst part, we wish to apply the electronic stress tensor 
analysis to the other elements which are conventionally classified as metalloids, B, Si and 
As, to see whether it can provide a criterion to dehne metalloids. For that purpose, it is 
also necessary to extend the analysis to the elements nearby metalloids in the periodic table. 
We also wish to apply the electronic stress tensor analysis to transition metals, ionic bonds, 
hypervalency, and so on, to strengthen its usefulness. It would enable us to deepen our 
understanding of the nature of chemical bonding. As for the second part, a further direction 
of the study will be to provide more evidence for our results by using other types of molecules. 
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If the relation concerning the force constant, which is an experimental observable derived 
from a vibrational spectrum, is established in general, it would help to solidify such ideas as 
the stress-tensor-based energy density and Lagrange surface. In the end, we would like to 
emphasize that our studies are based on the quantities dehned at each point in space, which 
originate from the quantum held theoretic consideration and not from the electron density. 
We believe that our method will lead us to new and benehcial views on chemical systems. 
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TABLE I: Data for the bonds between Ge, Sb, and Te atoms in the GST molecules, rg is the bond 
distance, and are three eigenvalues of the electronic stress tensor density at the 

Lagrange point. 


Molecule 

Te [A' 

.^S'33(x103) 

rg^22(xio3) 


GeHs-GeHs 

2.343 

-7.533 

-18.29 

-20.08 

GeH 2 (tBu)-GeH 2 (tBu) 

2.489 

-1.811 

-16.21 

-16.31 

GeH3-GeH2 

2.489 

-1.768 

-15.71 

-15.95 

GeH3-GeH2(tBu) 

2.486 

-1.841 

-16.23 

-16.28 

GeH3-GeH3 

2.483 

-1.861 

-16.25 

-16.25 

SbH2-SbH2 

2.948 

0.136 

-9.121 

-9.562 

SbH(iPr)-SbH 2 

2.946 

0.053 

-9.204 

-9.654 

SbH(iPr)-SbH(iPr) 

2.937 

-0.061 

-9.392 

-9.891 

Sb(iPr) 2 -SbH 2 

2.944 

0.022 

-9.335 

-9.724 

Sb(iPr) 2 -SbH(iPr) 

2.939 

-0.055 

-9.451 

-9.897 

Sb(iPr) 2 -Sb(iPr )2 

2.954 

0.213 

-9.226 

-9.752 

TeH-TeH 

2.854 

1.897 

-12.43 

-12.45 

Te(iPr)-TeH 

2.844 

1.514 

-12.69 

-12.82 

Te(iPr)-Te(iPr) 

2.831 

1.139 

-13.03 

-13.24 

GeH2-SbH2 

2.703 

-0.940 

-12.04 

-12.61 

GeH 2 -SbH(iPr) 

2.696 

-1.102 

-12.15 

-12.83 

GeH2-Sb(iPr)2 

2.692 

-1.218 

-12.25 

-12.99 

GeH2(tBu)-SbH2 

2.703 

-0.860 

-12.49 

-12.86 

GeH 2 (tBu)-SbH(iPr) 

2.702 

-0.894 

-12.59 

-12.97 

GeH2(tBu)-Sb(iPr)2 

2.704 

-0.869 

-12.58 

-13.01 

GeH3-SbH2 

2.700 

-0.811 

-12.52 

-12.82 

GeH 3 -SbH(iPr) 

2.694 

-0.995 

-12.70 

-13.04 

GeH3-Sb(iPr)2 

2.696 

-0.987 

-12.69 

-13.07 

GeH 2 (tBu)-Te(iPr) 

2.654 

0.059 

-14.36 

-15.28 

GeH 2 (tBu)-TeH 

2.658 

0.134 

-14.05 

-15.15 

GeH 2 -TeH 

2.653 

-0.649 

-13.82 

-15.08 

GeH 2 -Te(iPr) 

2.637 

-1.078 

-14.32 

-15.57 

GeH 3 -TeH 

2.649 

0.209 

-14.28 

-15.30 

GeH 3 -Te(iPr) 

2.636 

-0.183 

-14.76 

-15.75 

SbH 2 -TeH 

2.892 

0.838 

-10.60 

-11.12 

SbH 2 -Te(iPr) 

2.874 

0.343 

-11.05 

-11.52 

SbH(iPr)-TeH 

2.890 

0.627 

-10.59 

-11.28 

SbH(iPr)-Te(iPr) 

2.878 

0.441 

-10.97 

-11.58 

Sb(iPr) 2 -TeH 

2.892 

0.410 

-10.66 

-11.20 

Sb(iPr) 2 -Te(iPr) 

2.879 

0.338 

-11.05 

-11.53 
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TABLE II: Data for the bonds in the crystal structure of Ge, Sb, and Te. r^, is the distance between 
the nearest neighborhood atoms, and are three eigenvalues of the electronic stress 

tensor density at the Lagrange point. 


Crystal Ve [A] 

(xl03) rf2 

(xlO^) 

(xlO^) 

Ge 2.450 

-2.612 

-19.52 

-19.52 

Sb 2.907 

-1.981 

-12.12 

-12.23 

Te 2.835 

-1.712 

-17.46 

-17.81 
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TABLE III: Data for the bonds between Ge, Sb, and Te atoms in the GST molecules. B and k 
are respectively the bond dissociation energy and force constant, and our bond orders 

(Eqs. ()15ll and (|17l) i. a is a fitting parameter for the energy density distribution on the Lagrange 
surface (see the text for the details). 


Molecule 


B [kcal/mol] k [a.u.] ^e{s) 


GeHa-GeHa 

40.40 

0.122 

0.179 

0.648 

1.561 

GeHa(tBu)-GeH 2 (tBu) 

59.79 

0.080 

0.132 

0.485 

1.443 

GeHg-GeHa 

40.30 

0.079 

0.128 

0.473 

1.444 

GeH 3 -GeHa(tBu) 

59.86 

0.082 

0.132 

0.484 

1.448 

GeHg-GeHg 

59.58 

0.083 

0.132 

0.482 

1.455 

SbHa-SbHa 

32.21 

0.058 

0.071 

0.284 

1.186 

SbH(iPr)-SbHa 

32.62 

0.056 

0.072 

0.291 

1.175 

SbH(iPr)-SbH(iPr) 

33.09 

0.060 

0.074 

0.301 

1.177 

Sb(iPr)a-SbH 2 

32.72 

0.056 

0.073 

0.297 

1.166 

Sb(iPr) 2 -SbH(iPr) 

33.19 

0.056 

0.075 

0.303 

1.174 

Sb(iPr) 2 -Sb(iPr )2 

31.05 

0.056 

0.072 

0.294 

1.164 

TeH-TeH 

33.24 

0.064 

0.088 

0.288 

1.279 

Te(iPr)-TeH 

34.86 

0.063 

0.092 

0.307 

1.278 

Te(iPr)-Te(iPr) 

35.32 

0.066 

0.097 

0.327 

1.285 

GeHa-SbHa 

28.06 

0.069 

0.098 

0.378 

1.323 

GeH 2 -SbH(iPr) 

28.94 

0.068 

0.100 

0.390 

1.316 

GeH 2 -Sb(iPr )2 

30.10 

0.069 

0.102 

0.396 

1.317 

GeH2(tBu)-SbH2 

45.95 

0.060 

0.101 

0.389 

1.314 

GeH 2 (tBu)-SbH(iPr) 

45.75 

0.060 

0.102 

0.395 

1.309 

GeH 2 (tBu)-Sb(iPr )2 

46.22 

0.062 

0.102 

0.398 

1.298 

GeHg-SbHa 

45.13 

0.063 

0.100 

0.385 

1.323 

GeH 3 -SbH(iPr) 

45.72 

0.064 

0.103 

0.399 

1.318 

GeH3-Sb(iPr)2 

45.95 

0.065 

0.103 

0.401 

1.309 

GeH 2 (tBu)-Te(iPr) 

49.93 

0.064 

0.114 

0.404 

1.357 

GeH 2 (tBu)-TeH 

52.67 

0.064 

0.112 

0.415 

1.355 

GeHa-TeH 

33.96 

0.075 

0.114 

0.403 

1.393 

GeH 2 -Te(iPr) 

33.70 

0.077 

0.119 

0.428 

1.391 

GeHg-TeH 

50.13 

0.069 

0.113 

0.404 

1.375 

GeH 3 -Te(iPr) 

49.04 

0.071 

0.118 

0.427 

1.375 

SbHa-TeH 

35.71 

0.061 

0.080 

0.289 

1.246 

SbH 2 -Te(iPr) 

34.54 

0.052 

0.085 

0.310 

1.260 

SbH(iPr)-TeH 

38.52 

0.061 

0.082 

0.297 

1.236 

SbH(iPr)-Te(iPr) 

37.04 

0.061 

0.085 

0.297 

1.233 

Sb(iPr) 2 -TeH 

39.79 

0.048 

0.082 

0.301 

1.239 

Sb(iPr) 2 -Te(iPr) 

38.36 

0.060 

0.085 

0.319 

1.233 
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TABLE IV: Data for the homonuclear diatomic molecules, r^, oj, and B are respectively the bond 
distance, frequency, and dissociation energy, whose values are taken from the references. The force 
constant k is converted from u. Details are found in the text, and &£( 5 ) are our bond orders 
(Eqs. ([I5]) and (fTTll L 


Molecule 


[A] 


u cm 


-1 


fc [a.u.] i?[kcal/mol] References 


MS) 




0.74144 4401.21 0.367 


103.26 


62 


1.000 1.000 


Li2 

2.6729 

351.43 

0.016 

24.2 

162] 

0.013 

0.112 

Be2 

2.453 

270.7 

0.012 

2.308 


0.035 

0.165 

B2 

1.590 

1051.3 

0.230 

69 

[62] 

0.452 

1.205 

C2 

1.2425 

1854.71 

0.781 

143 

[62J 

2.349 

3.542 

N2 

1.09769 

2358.57 

1.474 

225.0 

[£2] 

5.849 

6.437 

O2 

1.20752 

1580.19 

0.756 

118.0 

m 

3.229 

4.079 

F2 

1.41193 

916.64 

0.302 

36.94 

\m 

1.303 

1.425 

Na2 

3.0789 

159.125 

0.011 

16.6 

m 

0.006 

0.065 

Mg2 

3.891 

51.12 

0.001 

1.16 

m 

0.004 

0.022 

A 12 

2.466 

350.01 

0.063 

37 

[62] 

0.074 

0.374 

Si 2 

2.246 

510.98 

0.138 

74.0 

[£2] 

0.245 

0.802 

P2 

1.8934 

780.77 

0.358 

116.1 


0.650 

1.701 

S2 

1.8892 

725.65 

0.319 

100.76 

[62] 

0.701 

1.513 

C 12 

1.988 

559.78 

0.208 

57.1742 

[6^ 

0.421 

0.645 

K2 

3.9051 

92.021 

0.006 

11.9 

[62] 

0.004 

0.044 

Ca2 

4.2773 

64.93 

0.003 

3.0 

m 

0.004 

0.038 

Ga2 

2.746 

162 

0.034 

26.36 

[65, 6^ 

0.042 

0.133 

Ge2 

2.3680 

287.9 

0.116 

65 

[62, £2 

0.122 

0.478 

As2 

2.1026 

429.55 

0.262 

91.3 

[£2] 

0.423 

1.286 

Se2 

2.166 

385.303 

0.225 

78.66 

[62] 

0.379 

1.015 

Bra 

2.2811 

325.321 

0.158 

45.444 

[62J 

0.246 

0.540 

Rba 

4.2099 

57.781 

0.005 

11 


0.004 

0.042 

Sra 

4.67174 40.32831 

0.003 

3.036 

[£21 

0.003 

0.031 

Ina 

3.14 

111 

0.027 

17.8 

[TO, 7^ 

0.034 

0.170 

Sna 

2.746 

189.74 

0.082 

43.8300 

[72, 73j 

0.108 

0.499 

Sba 

2.476 

269.623 

0.167 

69.0672 

[74] 

0.220 

0.869 

Tea 

2.5574 

247.07 

0.150 

61.73 

[62] 

0.195 

0.678 

I 2 

2.666 

214.50 

0.111 

35.5672 

[£2] 

0.137 

0.385 
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(c) Te(iPr )2 



(a) GeH 3 (tBu) 


FIG. 1: The optimized structures of the precursors: (a) GeH 3 (tBu) (Cat,), (b) Sb(iPr )3 (Ca), and 
(c) Te(iPr )2 (C 2 ). 
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FIG. 2: The largest eigenvalue of the electronic stress tensor density (color map) and corresponding 
eigenvector (black rods) are shown in panel (a), and tension density (black arrows) and its norm 
(color map) are shown in panel (b), around a G-G bond in GeH 3 (tBu). In both panels, the electronic 
interface is depicted by the green dashed lines, and the Lagrange point is shown by a black square. 
In panel (a), the red solid lines show the zero contour lines of the eigenvalue. 
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FIG. 3; The largest eigenvalue of the electronic stress tensor density and corresponding eigenvector 
between the GST atoms are shown in the same manner as Fig. [2] (a), (a) GeHs-GeHs, (b) SbH 2 - 
SbH 2 , (c) TeH-TeH, (d) GeH 3 -SbH 2 , (e) GeHa-TeH, and (f) SbH 2 -TeH. 
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FIG. 4: The relation between bond distance and the largest eigenvalue of the electronic stress tensor 
at the Lagrange point. The labels with “bulk” in the subscripts denote that they are computed for 
the crystal structures. 
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FIG. 5: Differential eigenvalues of the electronic stress tensor at the Lagrange point. In the panel 
(a), labels with “bulk” in the subscripts denote that they are computed for the crystal structures. 
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FIG. 6: The relation between energy density based bond order and (a) bond distance, (b) 
dissociation energy, and (c) force constant, for the GST bonds. 
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FIG. 7: In panel (a), the tension density is plotted similarly to Fig. [2] (b) for the Ge-Sb bond in 
a GeH 3 -SbH 2 molecule. The Lagrange surfaces are shown by black solid lines and the Lagrange 
points by black squares. In panel (b), the energy density distribution on the Lagrange surface 
between the Ge and Sb atoms is plotted as a function of the distance from the Lagrange point. 
The energy density is normalized by the energy density of a H 2 molecule at its Lagrange point. 
The black dashed line is the best fit Gaussian function: 6;; exp Q;|r — where a = 1.32. 
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FIG. 8: The ratio of b^(s) to for the GST bonds. 
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FIG. 9: The relation between force constant k and (a) and (b) ^e( 5 ), for the bonds in the 
homonuclear diatomic molecules. 
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